Differential light scattering: probing the sonoluminescence collapse 
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We have developed a light scattering technique based on differential measurement and polariza- 
tion (differential light scattering, DLS) capable in principle of retrieving timing information with 
picosecond resolution without the need for fast electronics. DLS was applied to sonoluminescence, 
duplicating known results (sharp turnaround, self-similar collapse); the resolution was limited by 
intensity noise to about 0.5 ns. Preliminary evidence indicates a smooth turnaround on a < 0.5-ns 
time scale, and suggests the existence of subnanosecond features within a few nanoseconds of the 
turnaround. 

PACS numers: 78.60.Mq, 42.65.Re, 42.68. Mj, 43.25.+y 
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Since Gaitan's seminal work [EQ, significant ad- 
vances have been made in our understanding of single- 
bubble sonoluminescence (SL). The dissociation hypoth- 
esis (DH) introduced by Lohse et al. j|| combines the 
merits of an intuitive approach based on the relatively 
tractable Raylcigh-Plcsset equation (RPE) of bubble dy- 
namics with an impressive ability to reproduce a wide 
range of observables [Qj5|^,^,||. More sophisticated the- 
ories f|,[l(| yield a more realistic picture of the phe- 
nomenon that is in general agreement with the results 
of the DH-RPE treatment. Experimentally, however, in- 
formation on the bubble interior is still very scarce. In 
particular, no direct and conclusive evidence exists yet 
for either plasma formation or shock waves inside the 
collapsing bubble. As a result, many competing theo- 
ries 0,[|jl|JT|,|l|,|l|,0|8| still vie with the adiabatic 
or shock-wave heating theory for the distinction of ac- 
curately describing the SL phenomenon. It seems de- 
sirable, then, to explore additional ways of probing the 
interior of the bubble with a time resolution compara- 
ble Jl9| to the duration of the flash, measured to be 40- 
380 ps [pO |2l|j2^ , |23t . Light scattering has already been 
shown |24[ 25| , |2t| to be a useful probe of the bubble dy- 
namics, sensitive as it is to the dielectric interface at the 
bubble wall. It is also a promising candidate for detec- 
tion of either plasma or shock waves, since both features 
can modulate substantially the local dielectric constant. 

Our goal was to push measurements of the light scat- 
tering cross section of the collapsing bubble to a greater 
temporal resolution than that afforded by the pulsed Mie 
scattering technique [2q| , which appears to be limited to 
around ^ ns by low light levels and the need for averaging. 
In order to achieve higher resolution, we have developed 
a technique called differential light scattering (DLS) that 
reduces statistical uncertainty in the detection process by 
making use of more powerful ultrafast laser pulses. Since 
DLS does not rely on a fiducial time reference, it, too, is 
completely insensitive to electronic timing noise, allowing 
the use of relatively slow detectors. The DLS technique 
is based on two central concepts: (i) using a differential 
measurement to yield jitter-free timing information, and 




FIG. 1. DLS concept applied to a simple curve, (a) Sam- 
pling of I(t) occurs in pairs; (b) for each pair, the difference 
signal is plotted against the hrst signal. Equal symbols in- 
dicate the correspondence between (a) and (b). Notice the 
slanted "wall" in (b) resulting from the two samples lying on 
opposite sides of the sharp cusp in (a). 

(ii) using polarized light to generate such a measurement 
through scattering. 

The differential measurement concept was recently in- 
troduced for the first time by Rella et al. in the con- 
text of ultrafast gating of optical pulses ^7|. The tech- 
nique they invented, differential optical gating (DOG), 
was used to measure the shape of a midinfrared pulse 
with subpicosecond resolution. Our technique applies the 
DOG concept to light scattering. DLS relies on collect- 
ing many pair of correlated samples of the same periodic 
event I(t) (see Fig. |l|). Each pair i consists of the first 
sample I(ti) and the second sample I(U + St), where ti 
is the time of the first sample (modulo the period T), 
and 8t is an appropriately chosen (and short) time delay. 
From each such pair, an intensity difference 



Sli = I{U + St) - I(U) 



(1) 



is produced and plotted against the first sample I(ti), 
generating what we will call a DLS plot. When enough 
event pairs are collected, the points representing them in 
the DLS plot will join together in defining a continuous 
curve. 
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FIG. 2. Experimental setup. Thin lines, electrical signals; 
thick lines, beam paths; dashed line, beam going into the 
page. 



The DLS plot can be thought of as a predictor: given 
an intensity / at some time t, the plot yields what / will 
be after a time St. One can numerically step along the 
curve on this plot to retrieve the desired direct function 
I(t). Depending on the nature of the features of interest, 
in some cases it is more fruitful to plot, for each pair, 8h 
against the second sample I(ti + St). The reconstruction 
of I(t) is then carried out backwards in time. When the 
data are particularly noisy, however, features apparent in 
the DLS plot will be lost in the reconstruction process, 
eliminating any benefit of the technique. In this case, it is 
better to work directly in DLS space. Using a model for 
I(t), a DLS curve can be generated from it by applying 
map ([!]), then fit to the data points in the DLS plot with 
a minimization algorithm. 

To implement the DLS concept (Fig. ||), a laser pulse 
is split equally in two and recombined as in a Michelson 
interferometer, with one pulse having traveled a longer 
path. The two-pulse train is focused onto the bubble, 
yielding two bursts of scattered light separated by an ad- 
justable time delay St. In order to distinguish between 
the two pulses during detection, an additional degree of 
freedom is needed. Color discrimination is a possibil- 
ity, but with several drawbacks, among which the need 
for frequency doubling and the strong wavelength depen- 
dence of light scattering. Polarization, on the other hand, 
is perfectly suited to this technique. Calculations using 
Mie scattering theory |p8| , which is rigorously valid for 
spheres of arbitrary size, show that scattered intensities 
are highly polarization dependent. For the case of lin- 
early polarized light, scattering at 9 — 90° (where zero 
is forward) vanishes if the polarization vector e is par- 
allel to the scattering plane. Therefore, the polarization 
of one of the pulses is made to rotate by 90° (with the 
quarter-wave plate shown in Fig. ||) before the two pulses 
are recombined. Figure ^ shows the sequence of the two 
pulses scattering from the bubble. The first pulse scat- 
ters preferentially in the plane containing one photomul- 
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FIG. 3. Scattering and polarization geometry. Not shown: 
lenses focus the input beam onto the bubble and relay scat- 
tered light onto the detectors. 



tiplier tube (PMT), while the second pulse does so in the 
perpendicular plane, which contains the other PMT. 

Combining the time delay and the polarization de- 
pendence results in the ability to assign the scattered 
intensity recorded by one detector to the earlier pulse, 
and that recorded by the other detector to the later 
pulse. Therefore, each laser burst yields an ordered pair 
of scattered intensities. By scanning the arrival time of 
the pulse pairs over some portion of the acoustic cycle, 
enough data can be collected to generate a DLS plot of 
the desired interval of the bubble's evolution. 

The experiments were performed in a 100-ml spheri- 
cal boiling flask filled with distilled, deionized water (re- 
sistivity p = 15.9 Mft-cm) at 20 ± 1 °C. The water was 
prepared in a gas-handling system under an air pressure 
of 0.2 bar, then loaded into the flask without further ex- 
posure to air. A sealed connection to a volume reservoir 
(kept at atmospheric pressure for these experiments) in- 
spired by Ref . p9[ provided pressure release from volume 
changes induced by temperature fluctuations. The entire 
assembly was leak tested; with the flask under vacuum, 
the rate of pressure rise was conservatively determined to 
be 16 nbar s . However, under normal operating con- 
ditions the flask is filled with water and repressurized to 
1 bar: this forces outside air to diffuse into the undersat- 
urated water through microscopic interfaces, resulting in 
a substantially lower rate of contamination. The experi- 
ments described here took place 110 days after loading. 
Within that time, the pressure in an initially empty flask 
would have risen to roughly 0.2 bar; the air concentration 
in the water-filled flask can instead be expected to have 
risen by perhaps a few percent of atmospheric saturation. 

The first acoustic resonance of the flask was deter- 
mined to be 26.9 kHz. The acoustic drive was provided 
by an audio amplifier, the output of which (typically 
~ 4.5 W rms) was fed through an impedance-matching 
network before being delivered in parallel to two disc- 
shaped piezoelectric transducers (PZTs) epoxied to dia- 
metrically opposite points on the flask. A third, smaller 
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PZT cemented to the flask provided acoustic pickup, used 
to map the normal modes of the flask and to monitor the 
behavior of the bubble through its filtered acoustic sig- 
nature. 

The laser used to probe the bubble was a regener- 
ative Tksapphire amplifier pumped by a Q-switched, 
frequency-doubled neodymium-doped yttrium lithium 
fluoride (Nd:YLF) laser operated at 1 kHz, 10 mJ/pulse 
(Positive Light Spitfire and Merlin, respectively) and 
seeded by a 82-MHz mode-locked TkSapphire oscillator, 
in turn pumped by an Ar + cw laser (Spectra Physics 
Tsunami and Bcamlok 2080, respectively). The oscil- 
lator provided 60-fs, 800-nm pulses at 82 MHz; the 
amplifier output consisted of partially uncompressed 
(chirped) 50-ps, 800-nm pulses at 1 kHz with approxi- 
mately 1 mJ/pulse. The dominantly TEMoo mode beam 
was sent through a spatial filter to clean up mode asym- 
metries and yielded a nearly Gaussian profile. To elimi- 
nate gross beam distortion caused by the irregular flask 
surfaces, a laser beam input port was made by cutting 
a hole in the flask and cementing in place a custom- 
made fused silica powerless meniscus. The light scattered 
by the bubble was collected with a relay system, passed 
through polarizers (appropriate for each branch) and 800- 
nm narrow bandpass filters, and delivered to two PMTs 
(Hamamatsu R955P and R636). The PMT signals were 
integrated by SRS SR250 boxcar averagers, which were 
in turn sampled by a 1-MHz A/D board on a personal 
computer. 

The synchronization scheme (shown schematically 
in Fig. |2|) involved generating a logic signal at 
facous = 26.9 kHz, and digitally dividing its frequency by 
27 to yield another logic signal at approximately 1 kHz. 
The 26.9 kHz logic signal was filtered before being fed to 
the audio amplifier to serve as the acoustic drive, while 
the 1 kHz signal was used to trigger the laser and the data 
acquisition electronics. This ensured that the SL drive 
signal and the regenerative amplifier pulse trains would 
be synchronized to each other to about 1-ns precision. 
Additional timing circuitry allowed for the delay between 
the SL flash (which occurs very nearly at the same point 
of the acoustic cycle, within 0.5 ns of turnaround |25|) 
and the laser pulse pairs to be varied continuously by up 
to 50 /is, either manually or automatically. This allowed 
us to probe the bubble at any given phase of the acoustic 
cycle. 

In order to obtain values for the ambient bubble radius 
Rq and the acoustic drive amplitude P a , we developed a 
time-stamp technique that yielded a time series of scat- 
tered intensities I(t) over the whole acoustic cycle. A 
time-to-amplitude converter (TAC, 566 EG&G ORTEC) 
measured the interval (up to a constant offset) elapsed 
between the arrival of the laser pulse and the SL flash, as 
signaled by an additional PMT sensitive to SL light only. 
The TAC output was logged through a boxcar along with 
the signal from one of the PMTs used in DLS, and used 
for time-stamping. Scattering events were recorded as 
the delay between the laser and SL was scanned auto- 
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FIG. 4. DLS plots of background-corrected data from a col- 
lapsing SL bubble. Data in (a) is uniformly scaled to account 
for higher attenuation than in (b). Overall range spanned by 
the data: (a) 90 ns; (b) 30 ns. The horizontal section A in 
(a) shows self-similar collapse behavior [ p6|plUl9[ ] lasting to 
within 5 ns of turnaround. The junction between sections A 
and B appears sharp in (a) but curved in (b), suggesting a 
smooth turnaround on a < 0.5 ns time scale. A barely re- 
solved kink in section B of (a) might be due to a subnanosec- 
ond feature. 



matically through a whole acoustic cycle. 

Since the result of this procedure was a time series 
of intensities, a calibration was performed to establish 
a conversion from I(t) to R(t). This was done using a 
stroboscopic imaging system similar to that of Ref . [p0[ , 
except that in our case the drive for the LED was locked 
to the same frequency / mus as that driving the bub- 
ble. We obtained R ma x by fixing the LED time delay 
so that the bubble was shown on the monitor screen at 
maximum size, and I m ax from the scattering data. A 
calculation based on Mie theory [^8| provided the I(R) 
map necessary to complete the calibration. 

In practice, uncertainties in the calibration of the imag- 
ing system, as well as in the actual measurement of the 
bubble size, prompted us to use our measurement of 
Rmax as an estimate with ± 10% uncertainty. The R(t) 
data were then fed to a fitting algorithm that established 
i?o, P a , and an appropriate overall scale factor in a non- 
linear least-squares calculation using the RPE. The fact 
that the scale factor for the best fit was determined in this 
way to be 1.09 ± 0.05 gave us confidence in the validity of 
our imaging method. The bubble parameter values thus 
found were Rq = 5.3 ± 0.2 fjjn and P a = 1.34 ± 0.04 bar. 

It is worth mentioning that the time-stamp technique 
described above can be used directly to obtain light scat- 
tering data from a collapsing SL bubble. The drawback 
is that, unlike in DLS, the electronic response of the mea- 
suring instruments is the limiting factor. We used this 
procedure to collect rough timing information as a cross- 
check in our analysis of DLS data; with the devices at 
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our disposal, 2-ns resolution was achieved. We estimate 
that with two microchannel plate PMTs, two constant- 
fraction discriminators, and a faster TAC, an overall tim- 
ing uncertainty of 50 ps should be achievable EO) . 

In Fig. ^ we show representative results from our DLS 
experiments. In these plots, as in Fig. [j] (b), the abscissa 
is I(U) and the ordinate is 5Ii. In Fig. |j(a) the delay 
between pulses was 5 ns, and in Fig. |(b) it was 1 ns. The 
range of 5t for which useful information can be gathered 
is dictated by the physical process under study: delays 
much shorter than 1 ns yielded DLS plots unresolved into 
a discernible structure, while delays much longer than 
5 ns are not well suited for investigating short time scales. 

To aid interpretation, we divide the plots into three re- 
gions: A (the collapse, t < 0), B (the transition region), 
and C (the rebound, t > 0). The approximately flat 
region A corresponds to the collapse, since 8Ii < 0, indi- 
cating that the bubble is shrinking. In C the rebounding 
bubble expands, but at a much lower rate than during 
collapse, so 8Ii is positive and smaller in magnitude than 
in A. However, a greater spread in the data there results 
in fuzzy clustering across Sli — 0. The straight "wall" 
in region B forms when the two pulses straddle t = 
(compare to the open squares in Fig. FP). 

The straight section A in Fig. |I](a) is due to a con- 
stant slope of I(t) during collapse. This critical behav- 
ior 1 31 19] has been previously observed for time scales 
ranging from 1 /is to 20 ns prior to turnaround p6| ; our 
observations extend it to t = —5 ns. In Fig. [|(a) sections 
A and B join rather abruptly, indicating a sharp cusp in 
I(t) on the time scale of the measurement (5 ns). This 
was expected given the measurements in Ref. . In Fig. 
^(b), however, section A appears to show a slight upturn 
before joining section _B, indicating a smooth transition 
on a time scale less than the pulse delay of 1 ns (since 
the "wall" section B is still discernible). Scatter in the 
data prevents a conclusive interpretation, but the avail- 
able evidence would support an estimate of the bubble 
turnaround time at a few hundred picoseconds. 

The collection lenses used have a /-number of 1.5; the 
finite acceptance cone they subtend introduces a pollu- 
tion, or cross talk, of unwanted light from the other pulse 
in each detector. Because of the strong polarization de- 
pendence of scattering, this cross talk is quite small: it 
was calculated from Mie theory, and confirmed experi- 
mentally, to be less than 5% of the total scattered inten- 
sity. Electrical cross talk was measured to be less than 
5%. The resulting overall intensity uncertainty is there- 
fore around 7%; the difference uncertainty varies across 
the plot. While in sections A and B the error estimates 
are consistent with the observed spread in the DLS data, 
in section C the spread is significantly larger. This has 
been observed before in Xe-filled bubbles and ascribed 
to nonsphericity ]2q| ; such asymmetry is reported here 
as regularly occurring in air-filled bubbles. 

In conclusion, we have introduced DLS, a light scat- 
tering technique based on the DOG |^7j concept of dif- 
ferential measurement and on sensitivity to polarization 



that uses intense ultrashort laser pulses to bypass the 
problem of electronic timing jitter. The intensity spread 
in the data is currently the limiting factor in the reso- 
lution achieved with this technique. Effectively, inten- 
sity noise is translated into timing noise by the mapping 
that a DLS plot generates. Accordingly, the resolution 
in the data shown is approximately 0.5 ns. The intrinsic 
resolution of DLS, however, is given by the laser pulse 
width used: with our equipment that can be made as 
low as 500 fs. Data collected from a collapsing SL bub- 
ble confirm earlier findings of a self-similar solution and 
of subnanosecond turnaround time; our preliminary re- 
sults suggest that the turnaround is smooth on a time 
scale of a few hundred picoseconds. 
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